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ABSTRACT 

In a series of recent papers, a new formalism has been developed that explains the inner structure of 
dark matter halos as collisionless, dissipationless systems assembled through mergers and accretion at 
the typical cosmological rate. Nearby ellipticals are also collisionless, dissipationless systems assem- 
bling their mass through mergers, but contrarily to the former structures they do not continuously 
accrete external matter because they are shielded by their host halos. Here we explore the idea that 
the infall of their own matter ejected within the halo on the occasion of a violent merger can play a 
role similar to external accretion in halos. The predicted stellar mass density profile fits the observed 
one, and the empirical total mass density profile is also recovered. 

Subject headings: gravitation — galaxies: formation — galaxies: structure — galaxies: elliptical and 
lenticular — galaxies: halos — dark matter 
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1. INTRODUCTION 

Nowadays there is increasing evidence about the im- 
portant role of major, non-dissipative mergers (Toomre 
& Toomre 1972) in the mass assembly of elliptical galax- 
ies (Conselice 2003; Bell et al. 2005) at intermediate or 
low-redshifts, the formation of their stellar populations 
and the dissipative processes previous to it having oc- 
cured mainly at high-z (Faber et al. 2005; de Lucia et 
al. 2005; Dominguez-Tenreiro et al. 2006). The prob- 
lem with dissipationless major mergers is that there is 
no known analytica l treatment for the violent relaxation 
(|Lvnden-Belllll967h they lead to. 

Elliptical galaxies are similar to cold dark matter 
(CDM) halos in many respects. Apart from the sim- 
ilar role played by major, non-dissipative mergers in 
their mass assembly, their mass, velocity dispersion, and 
length scales show systematic regularities and correla- 
tions that can be related with each other (Bernardi et 
al. 2003; Graham et al. 2006; Ohorbe et al. 2007). Like- 
wise, their 3D density profiles are all well fit (Navarro et 
al. 2004; Merritt et al. 2005; Merritt et al. 2006; Onorbe 
et al. 2007) by the Einasto (Einasto & Haud 1969) law 
or the Prugnicl-Simien (1997) approximate analytical in- 
version of the projected Sersic (1968) law. These co- 
incidences suggests that these profiles could have been 
shaped by similar physical processes. 

In a series of recent papers, a new formalism has been 
developed (Salvador-Sole et al. 2007, hereafter SMGH 
and references therein) that explains the inner structure 
of halos from the very collisionless, dissipationless na- 
ture of CDM, and which predicts halo structural and 
kinematic properties that are in good agreement with 
the results of A-body simulations (Gonzalez-Casado et 
al. 2007). Thus, it is natural to consider whether this 
formalism can also explain the inner structure of ellipti- 
cals. 
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But things are not that simple. Halos are permanently 
accreting external matter, while ellipticals are not. As 
numerical simulations show (e.g., Dominguez-Tenreiro et 
al. in preparation), matter falling into the halo is too 
energetic to stick to the center, and it is deposited at 
the edge of the system. Only clumps massive enough 
are braked by dynamical friction and spiral down until 
merging with the central galaxy. Therefore, apart from 
such discrete (minor and major) mergers there is no con- 
tinuous accretion into central ellipticals. 

This is a crucial difference in the context of the SMGH 
model because in that model the structure of halos is de- 
termined precisely by their continuous accretion. Indeed, 
the profile at the edge of a halo of any extensive property 
is set by the rates of infalling and rebounding matter, 
which do not depend on the particular inner mass distri- 
bution but only on the current accretion rate. As CDM 
is collisionless, the spatial distribution of any property 
is necessarily at all derivative orders, implying that the 
whole respective inner profile adapts to the external one 
defined by accretion. As CDM is in addition dissipation- 
less, the steady inner region remains unaltered during 
accretion, causing the profiles to grow from the inside 
out. This specific growth allows us to infer the profile of 
any extensive property from its typical rate of increase by 
accretion (see below for the case of the density profile) . 

Although local ellipticals do not accrete external mat- 
ter, for a short interval after a merger they do collect the 
(bound) matter ejected in that violent event. The aim of 
the present Letter is to explore the possibility that such 
an infall of ejecta can play a similar role in ellipticals as 
that plaied by cosmological accretion in dark halos. 



2. DYNAMICS OF EJECTA 

When two galaxies pass each other close enough to be 
tidally disrupted their content is ejected in all directions. 
As the disruption takes place when the merging galax- 
ies are orbiting around the center of mass, the typical 
radius R p of the ejection region is larger than the half- 
mass radii r e of the progenitors. The particle velocities in 
the disrupted system should be approximately normally 
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distributed, 
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with both the proportionality factor and the characteris- 
tic 3D velocity dispersion <t p independent of the particle 
mass (the acceleration undergone by particles does not 
depend on it). 

Particles ejected at larger radial velocities expand more 
rapidly than those ejected at smaller ones. Furthermore, 
for a given velocity, particles ejected from inner regions 
experience a smaller gravitational pull and expand more 
rapidly than those ejected from outer regions. Therefore, 
ejecta tend to segregate into concentric shells with out- 
wards increasing radial velocities. Most of these shells 
follow bound orbits. They reach the apocenter and fall 
back, cross the central region and rebound to a some- 
what smaller turn-around radius owing to the crossing 
with shells falling in for the first time, and so forth. The 
result of this chaotic motion is a central relaxed object 
into which bound ejecta continue to fall and rebound for 
some time. 

These dynamics greatly resemble those of the cosmo- 
logical evolution of density perturbations, with the dif- 
ference that in the present case, apart from the shell- 
crossing between infalling and rebounding layers, there 
is some additional shell-crossing owing to the finite size 
of the ejection region as there tends to be segregation 
of particles according to radial velocities. In the present 
Letter, we do not consider these finite size effects. We fo- 
cus on the density profile of the steady object that would 
emerge were all particles ejected from the same typical 
radius i? p , and the above mentioned segregation would 
be instantaneously and perfectly achieved. In these cir- 
cumstances, some known results for the collapse of spher- 
ical density perturbations hold. 

If rebounding shells were suppressed, there would be no 
shell-crossing (see below) and the stellar mass M(v) in- 
side the shell with positive initial radial velocity v would 
be constant and equal to the total mass of shells with 
smaller velocities. Taking into account equation {l}, af- 
ter integrating over the tangential velocity components, 
we obtain for a radial velocity v 

M(v) = Aj dwexpf-— 1 ' 



if 
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(2) 

Similarly, the constant energy of stars inside that shell 
(for the system truncated at its radius and the potential 
origin at infinity) is 

w 2 cr 2 GM{i) 

y + y 



E(v) =A dv 
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The integrals in equations @ and ([3]) extend over posi- 
tive radial velocities only because particles with negative 
initial radial velocity, —v, cross the system, and join the 
shell with v, the time required to do so being one of the 
finite size effects we ignore. By taking v equal to infinity 
in these equations, we arrive at 



2E p . GM, 



— ^-2$ h (i? p ) and A = J £ 

rjRp V 7r cr p 
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relating the velocity dispersion <r p and the normalization 
constant A, which is equal to the proportionality fac- 
tor in equation |T]) times the average particle mass 1 , to 
the typical ejection radius R p and the total stellar mass 
Mp, energy E p , and baryon mass fraction rj of ejecta. 
In equations and (|3|) , G is the gravitational constant 
and 

$ h (r) = -i-irG Po rln h P [3n h , (r/r h ) 1/nb 
x (^r(3n h ) + T(2n h ) {l - P [2n h , (r/r h ) 1/,lh ] }) (5) 

is the potential of the steady dark halo, endowed with an 
Einasto density profile (SMGH), 



p{r) = po exp 



l/nh 



(6) 



with two independent shape parameters, say, ny, and r n 
for one fixed total mass In equation (O, r(x) is the 
gamma function and P(a, x) is defined as 



P(a,x) 



r( 
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In deriving equation |3j), we have taken into account that, 
as the CDM distribution in the steady halo adapts to ex- 
ternal accretion, the void left by the ejected dark matter 
will be rapidly filled again, causing the net amount of 
ejected CDM to decrease. For this reason, rj may be sub- 
stantially larger than the original baryon mass fraction 
r]p within R pi the relation between the values being 



n 



n^Mp - M h (R p ) ' 
where M p is the total ejected stellar mass and 



(8) 



(9) 



is the steady halo mass inside r for the Einasto density 
profile (eq. [B]). 

As mentioned above, rebounding and shell-crossing 
lead to a steady state object, which extends, at any given 
time t, out to the limiting radius R(t). This radius can be 
estimated from the virial relation W(t) — 2E(t), where 
W(t) is the potential energy of the relaxed stellar system, 
truncated at R(t) and assumed to have a uniform mass 
distribution, and E{t) the corresponding conserved (in 
the previous model with no shell-crossing) total energy. 
Estimating E(t) is at turn-around and assuming a uni- 
form mass distribution, in the absence of any external 
gravitational potent ial, leads to R(t) e qual to half the 
turn-around radius (jGunn fc GottJ[l972T ). In the present 
case, E(t) is given directly by equation ([3]) where v is 
equal to the initial radial velocity of the shell collapsing 
at t. Dividing that virial relation by the corresponding 
inner stellar mass (eq. [2]), we are led to the following 
implicit equation for R(v) 



GM(v) 
rjR(v) 



2 2a 2 v 2GM(v) 



r\Rp 



2$ h (i? p ). 
(10) 



1 The particle mass factorizes because the velocity distribution 
function (TTJ does not depend on it. 
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Fig. 1. — Typical Prugniel-Simicn density profiles for ellipticals 
with stellar masses equal to 3 X 10 9 Mq (long-short-dashed line), 
10 10 M (dot-dashed line), 3 X 10 11 M Q (dashed line), and xlO 11 
Mq (dotted line) according to Graham et al. (2006 and references 
therein) and their best fit by the theoretical density profile derived 
here (solid lines). See Table 1 for the values of the respective 
parameters. 



The initial velocity v of the shell collapsing at t is sim- 
ply the inverse of the collapse time, in the model with 
no shell-crossing, of the shell that starts with v. This 
collapse time, t co u(v), can be obtained by numerical in- 
tegration of the equation of motion 



G^M^J + AfhtK*)]} 



(11) 



then imposing r(t) equal to zero. We have checked that 
t co \\{v) is an increasing function of v, which proves that 
there would be no shell-crossing if the rebounding shells 
were suppressed. 

The density profile of the steady object at t can be ob- 
tained in the way explained in SMGH provided the struc- 
ture is stable against the infall of ejecta which warrants 
its inside-out growth. In other words, the characteristic 
time of stellar mass collapse (in the model with no shell- 
crossing), inverse of the collapse rate lZ a = M/M, must 
be larger than the crossing time for particles at R(t) for 
every collapse time t or, equivalently, for every initial 
velocity v of bound shells, 



T]R 3 (v) 



-,1/2 



2GM(v) 



(12) 



where 



1 dM(v) /dtcoii 



M(v) dv \ dv 
Condition (TTZ|) is satisfied, so the SMGH model holds 
3. PREDICTED DENSITY PROFILE 
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Fig. 2. — Baryon mass fraction (top panel) and total (baryon or 
stellar plus CDM) density profiles (bottom panel) for the theoret- 
ical ellipticals (same symbols) plotted in Figure [l] 



Given the inside-out growth of the steady system, 
mass evolves according to the equation 



its 



M(t) 



R.(t) 



dr 4np(r) r 2 



(14) 



where the density p is independent of time. Differenti- 
ating equation ([l"4|) leads to the density profile p[R(t)\ 
in terms of the time derivatives of both the collapsing 
mass M(t) and the virial radius R(t). More simply, by 
expressing M (t) and R(t) as M[v(t)] and R[v(t)] and dif- 
ferentiating equation (fT4"|) by means of the chain rule, the 
time derivative of v simplifies out, and, after some alge- 
bra, we are led to the following parametric equations for 
the density profile we want, 



p[v) 



M(v) + r]M h [R(v)} 



4ttR 3 (v) 



1- 



2R(v) 



Ru 



7r r)R p a p v ( 3w 2 
1-J- exp — 



6 GM D 



(15) 



and R(v) given by equation (flQl) . 

In Figure [TJ we plot the typical 3D density profiles of 
the Prugniel-Simien form, 



p(r) 



p.. \ — ) exp { -d n 



l/n 

-I - 1 

r. 



. (16) 



for present day ellipticals of various stellar masses M. 
For p e equal to the density at the half-mass radius r e , d n 
and p become functions of n, approximately given, for 
n > 0.5, by 2n - 1/3 + 0.009876/n and 1.0 - 06097/n + 
0.05463/n 2 , respectively (see Merrit et al. 2006). Thus, 
the profile (|16[) is fixed by three independent parameters, 
reducing to two for one given mass M. The typical values 
of these parameters have been taken from Graham et 
al. (2006 and references therein); see Table 1. 
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(kpc) 


Mq pc" 3 ) 


(kpc) 
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0.3 


1.70 


0.94 


-0.94 


0.20 (0.05) 


6.7 (6.2) 


0.47 (0.09) 


2.4 (2.2) 


24.9 (21.7) 


60.0 (36.0) 


0.05 (0.16) 


1.0 


2.34 


1.12 


-0.70 


0.35 (0.13) 


7.0 (6.6) 


0.25 (0.71) 


2.5 (2.3) 


27.4 (18.6) 


92.0 (76.0) 


0.14 (0.27) 


3.0 


3.10 


1.57 


-0.72 


0.55 (0.42) 


7.2 (7.1) 


0.15 (0.20) 


2.6 (2.5) 


34.3 (34.7) 


126 (112) 


0.23 (0.28) 


10 


4.28 


3.16 


-1.17 


1.40 (1.40) 


7.7 (7.7) 


3.85 (3.85) 


2.9 (2.9) 


58.4 (58.4) 


197 (197) 


0.29 (0.29) 


30 


5.72 


10.0 


-2.26 


4.50 (4.25) 


8.4 (8.4) 


49.6 (55.4) 


3.3 (3.3) 


99.7 (108) 


355 (335) 


0.27 (0.28) 



TABLE 1 



For each empirical profile, we also plot in Figure [T] its 
best fit to the theoretical profile given by equations (115|) 
and (fTT)|) . As ejecta are rapidly collected by the newborn 
elliptical (most bound mass is collected in less than one 
Gyr), its typical stellar mass M has been assumed to 
coincide with the final asymptotic value, 



M = M p — A 



d,exp(- — 



M p erf 



V3 / GM P 
cr p V riR p 



$h(i? P ) 



1/2 



(17) 



where v csc = y/2GAl p / (r)R p ) — 2$h(i? p ) is the escape ve- 
locity at R p . Then, the typical total dark halo mass Mh 
has been obtained from M by means of the relation by 
Shanks et al. (2006). According to these authors, the ra- 
tio of dark halo to stellar mass, equal to ~ 14 for normal 
ellipticals (see also Ohorbe et al. 2007), rapidly increases 
towards the dwarf mass end. But the estimate is quite 
uncertain there. For this reason, we have also considered 
the assumption of an ever constant dark halo to stellar 
mass ratio equal to 14. Then, the Einasto shape pa- 
rameters fixing the CDM halo density profile have been 
obtained using the SMGH prescription ignoring any adi- 
abatic contraction. Having fixed all these values (see 
Table 1) there are only three free parameters to be ad- 
justed: <7 P , Rp, and r/ p . 

As can be seen from Figure [T] the theoretical profile 
yields a very good fit for very reasonable values of those 
free parameters (see Table 1): a p takes values of the order 
of (or slightly larger than) the typical velocity dispersion 
of the progenitors, R p takes values of the order of their 
limiting radii, and r/p of the order of their total baryon 



mass fraction, (see Fig. 13 in Ohorbe et al. 2007). 

Figure [2] shows the predicted baryon mass fraction and 
total (stellar plus CDM) density profiles, which also agree 
very well with empirical data (Koopmans et al. 2006; 
Ohorbe et al. 2007). Note, in particular, the isothermal- 
like behavior of the total density profiles. 

4. CONCLUSIONS 

These results seem to confirm that the structure of 
nearby ellipticals (and of simulated non-accreting CDM 
halos; see Hansen et al. 2005) is set in major mergers 
(minor mergers let it essentially unaltered), in the way 
explained in SMGH through the infall of matter ejected 
at that violent event. At this stage, we found no evidence 
of a different formation mechanism for dwarf ellipticals 
compared to normal bright ones. 

In the present Letter, all ejecta were assumed to be 
thrown from the same typical radius. This limits severely 
the minimum radius down to which the predicted den- 
sity profile can be calculated to only ~ R p /2 (exactly the 
value in the absence of a dark halo; see eq. [10]), corre- 
sponding to the virial radius for shells initially at rest (at 
turn-around). In a forthcoming paper, we will apply a 
more accurate treatment that will allow us to reach the 
galaxy center. In that paper, we will aso deal with the 
kinematics of ellipticals. 
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